Variation after parity projection calculation with Skyrme interaction for light nuclei 
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A self-consistent calculation with variation after parity projection is proposed to study both 
ground and excited states of light nuclei. This procedure provides description of the ground state 
incorporating some correlation effects, and self-consistent solutions for the excited states of nega- 
tive parity. For flexible description of nuclear shapes, single particle orbitals are represented on a 
uniform grid in the three-dimensional Cartesian coordinates. The angular momentum projection 
is performed after variation to calculate rotational spectra. To demonstrate the usefulness of the 
method, results are shown for two N = Z nuclei, 20 Ne and 12 C, for which clustering correlations 
are known to be important. In the 20 Ne nucleus, both cluster-like and shell-model-like states are 
described simultaneously in the present framework. For 12 C nucleus, the appearance of three-alpha 
clustering correlation in the ground state is investigated in relation to the strength of the two-body 
spin-orbit interaction. 



PACS numbers: 21.10.Hw, 21.60.Jz, 27.20.+n, 27.30.+t 



Nuclear mean-field calculations have been successful 
for systematic description of nuclear ground state prop- 
erties with a few adjustable parameters [lj. The mean- 
field theory has also been extensively applied to descrip- 
tion of excited states. For example, collective vibrational 
as well as single particle-hole excitations are described 
by the random phase approximation which is equivalent 
to a small amplitude limit of the time-dependent mean- 
field dynamics. Excited states may also be described as 
local-minimum solutions if there is more than one self- 
consistent solution in the static calculation. In principle, 
the generator coordinate method is expected to provide a 
unified description for any kind of excitations. However, 
in practice, one need to specify a few important generator 
coordinates in advance by physical intuition. 

In this article, we present an attempt to describe both 
ground and excited states in the self-consistent approach 
with the variation after parity projection (VAPP). For 
an even-even nucleus, the lowest energy positive-parity 
solution describes the ground state in which the correla- 
tion beyond the mean field is incorporated to a certain 
extent. Self-consistent negative-parity solutions should 
correspond to the negative-parity excited states. 

The theory of the variation after projection (VAP) of 
the symmetry-violating internal state has a long history 
0- However, practical applications with full variation 
of single-particle orbitals are rather few even for projec- 
tion with respect to the parity. Previously we achieved 
the VAPP calculation employing the uniform grid repre- 
sentation in the three-dimensional (3D) Cartesian coor- 
dinates in which the simple BKN interaction was used 
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In the ordinary mean- field calculations, one usually 
obtains self-consistent solutions with axial and reflection 
symmetries for most nuclei. However, in the VAPP calcu- 
lations, the self-consistent solutions are found to violate 
these symmetries. A preliminary report on the realistic 
calculations employing the Skyrme interaction has been 
presented in Refs. Qllg- 

We apply the method to light N = Z nuclei, 20 Ne 
and 12 C. There are several reasons why we consider the 
light nuclei. For light nuclei, the VAP may significantly 
modify the mean-field solutions and lead to the energy 
gain by the projection. In the excited states of light nu- 
clei, various cluster structures have been found to ap- 
pear. They have been successfully analyzed with the 
microscopic cluster models in which the VAP calcula- 
tion was carried out 6]. The present approach is useful 
to understand the mechanism for the appearance of the 
cluster structures. Finally, there are rapid progresses in 
the experimental research on light unstable nuclei. The 
present framework will be useful to analyze and predict 
structures of unstable nuclei. 

The antisymmetrized molecular dynamics (AMD) 
method developed by Kanada-En'yo and Horiuchi Q has 
been utilized to describe clustering phenomena in light 
stable and unstable nuclei. In the simplest version of 
AMD, the method is identical to the VAPP with a re- 
stricted Slater determinant in which each orbital is a 
Gaussian wave packet. This may be regarded as an ap- 
proximation to our present approach. The feasibility of 
the AMD allows them to perform the variation after an- 
gular momentum projection |8(- On the other hand, the 
present work brings more accurate description of the sin- 
gle particle orbitals, and establishes an intimate link to 
the Skyrme Hartree-Fock theory. 

We start with a brief description of the theory 
We consider an arbitrary Slater determinant $ = 



2 



■^^det{fa(xj)}, with x = (r,a), the space and spin co- 
ordinates. We apply the parity projection operator for 
this Slater determinant, = -^=(1 ± P)$, where P 

is the space inversion operator. We then consider the 
energy expectation value for the states with definite par- 
ity, and make a variation with respect to the single 
particle orbitals fa, 



usual, 
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Here we imposed two kinds of constraints in the above 
variation |3j|. The first one with the Lagrange multipli- 
ers eij is introduced to orthonormalize single particle or- 
bitals, fa. The second with the multiplier if is to coincide 
the center-of-mass of the wave function $ with the ori- 
gin of the coordinate, Y^ii'Pi\^\ ( t > i) = 0- This constraint 
assures that the parity operation is made with respect 
to the center-of-mass of the nucleus and minimizes the 
occurrence of the spurious center-of-mass excitations by 
the parity projection. 

The variation of Eq. yields the following self- 
consistent equation, 

{h - ff ■ r) fa ± {$\P\$}{h P fa - 

3 

+ (E^-E)4> i = J2^s (2) 



where h is the usual Hartree-Fock Hamiltonian. hp has 
the same structure as h, however, all the densities are re- 
placed with the transition densities which are the matrix 
elements of density operators between the wave function 
<!> and its parity- inverted state P<&. E^ is the energy 
expectation value with respect to the wave function 
while E is with respect to $. fa is defined by 



(3) 



with By = (fa\P\fa). 

In the practical calculations, Eq. J2J is n °t solved di- 
rectly. Instead, the imaginary-time method is employed 
in which the left-hand-side of Eq. J5J is used as the gradi- 
ent of the energy functional. We discretize the 3D Carte- 
sian coordinates into uniform square grid and represent 
the single particle wave functions on the grid points. The 
grid spacing is taken to be 0.8 fm, and the grid points in- 
side a sphere of radius 7.2 fm are used in the calculations 
below. 

After obtaining the self-consistent solutions, we make 
the angular momentum projection (AMP) to calculate 
rotational spectra. The self-consistent solutions in the 
VAPP are usually not axially symmetric. Therefore, one 
must perform three-dimensional rotation in Euler angles, 
fi = (a, (3, 7), for the AMP. We define the AMP state as 
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where is the rotation operator and D J MK (Q) is the 

Wigner's D-function. We then define relevant matrix 
elements, hi^g] and n^g}. The Hamiltonian matrix ele- 
ment h'xtf] is given by 
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KK' 
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4), (5) 



and the similar expression for the norm matrix element, 

71kk' ■ ^ or P rac tice, we achieved the rotation of the single 
particle orbitals by successive rotations of a small angle. 
Following the Taylor expansion method for time evolu- 
tion which was employed in the time-dependent Hartree- 
Fock calculations we calculate the rotation of small 
angle Aa around z-axis by 



e AQ = E 

k=0 



(Aa)* 
k\ 



(6) 



where N max is taken to be 4. Typically, each Euler angle 
is discretized into 200 steps for a and 7, and into 400 
steps for j3. In Eq. (JSJ , the integrations over Euler angles 
are achieved with discretization of 20 points for a and 
7 and 400 (all) points for j3. It is important to divide 
rotations of three Euler angles into two and one; the ro- 
tations of the angles j3 and 7 for the ket state, and the 
rotation of the angle a for the bra state, to reduce the 
computational costs. 

We employ the Skyrme interaction in constructing the 
energy functional in Eq. (Q. There is ambiguity in the 
choice of the density of density dependent force for off- 
diagonal matrix elements [T^j, which appear in Eqs. QJ, 
P)l. and (J5J. In this paper, we simply use corresponding 
transition densities for them. We will hereafter abbre- 
viate the present scheme as the parity-projected Skyrme 
Hartree-Fock (PPSHF) method. 

Before showing our results, we mention some prob- 
lematic aspects of the Skyrme interaction in the present 
framework which we have encountered in the practical 
calculations. Since the Skyrme interaction has been so 
constructed to be used for a single Slater determinant 
state, it is by no means evident whether it may provide 
a reasonable description in the VAPP calculations. In- 
deed, we have observed that, if we make a variation of 
single particle orbitals without any restriction, we ob- 
tain an unphysical solution in which the time-reversal 
symmetry is violated in the ground state. For exam- 
ple, the total binding energy of 20 Ne in the ordinary 
Skyrme Hartree-Fock (SHF) calculation with SGII force 
is about 170 MeV, while we obtain the binding energy of 
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TABLE I: Self-consistent solutions of Ne are summarized. 
The third row ("Energy (MF)") shows the results of ordinary 
SHF and PPSHF (IT = + , K* = 2~, K n = 0~) before 
AMP. Total energy is shown for the ground state (SHF and 
K n = + ) , while the excitation energies with respect to the 
ground state are shown for the excited states (K n = 2~ and 
K* = 0~). The fourth raw ("Energy (MF+AMP)") indicates 
the results for the lowest angular momentum state after the 
AMP. The fifth raw ("-Els") is the expectation value of the 
two-body spin-orbit interaction. The P2, 7, Pso and #32, are 
the density deformation of the internal single Slater determi- 
nant. 





K n = 0+ 


K* = 2~ 


K* = 0" 


SHF 


Energy (EXP) 


-160.645 


+4.97 


+5.78 


-160.645 


Energy (MF) 


-168.232 


+5.47 


+7.96 


-167.943 


Energy (MF+AMP) 


-171.333 


+4.91 


+6.42 




Els 


-8.517 


-15.350 


-5.443 


-8.775 




0.535 


0.589 


0.694 


0.728 


1 


0.0 


0.0 


0.0 


0.0 


/3so 


0.314 


0.0 


0.576 


0.0 


/3 32 


0.0 


0.178 


0.0 


0.0 


230 MeV with the 


PPSHF. 


Since the 


parameterization 



of the Skyrme interaction relevant to the time-odd com- 
ponent has not been fully tested, we restrict our wave 
function to be time-reversal invariant. Namely, we force 
nucleons either to fully occupy or unoccupy a pair of or- 
bitals which are mutually related to each other by the 
time-reversal operation in the Slater determinant. This 
restriction removes most part of the problems. However, 
in the positive parity solutions of some light nuclei, we 
still encounter a problem that the density of the self- 
consistent solution shows a small but unphysical oscilla- 
tion, namely, the density shows staggering in the neigh- 
boring grids. This problem can be overcome if we further 
ignore a part of the Skyrme energy functionals including 
terms p, j, T 0,0]. For a single Slater determinant 
with time-reversal invariance, these time-odd densities, 
p, j, T, vanishes identically. However, this is not the 
case for transition densities between <I> and P$. 

The problems are related to the fact that the superpo- 
sition of two Slater determinants by the parity projection 
may represent more varieties of correlation effects than 
the single Slater determinant may do. We leave it an 
interesting future problem to find appropriate Skyrme 
energy functional to be used in the PPSHF framework. 
In the following calculations, we will employ SGII param- 
eter set [l^. This force has been successful in describing 
the response properties of nuclei, although the absolute 
binding energies are not reproduced accurately. 

We first show the results of 20 Ne. In TableQ] the calcu- 
lated results are summarized. The definition of the defor- 
mation parameters is in Ref. 0|. Octupole deformation 
parameters foi and P33 are negligible and not shown. 
The positive-parity solution corresponds to the ground 
state. In the SHF calculation, this has a prolate shape 
with reflection symmetries. On the other hand, the PP- 



SHF calculation gives the solution with substantial P30 
deformation. The energy in the PPSHF calculation gets 
lower by 0.3 MeV from the ordinary SHF calculation. 

There are two solutions in the negative parity. The 
lowest energy solution has a small P32 deformation as 
well as the prolate deformation. The next solution has 
a large P30 deformation as well as /?2- In the angular 
momentum projection, we found that the former solution 
is characterized by = 2~ and the latter solution by 
K* = Q-. 

Previously we reported the existence of two local min- 
ima in the negative parity Q. After careful examina- 
tion, we have found that the K n = 0~ solution is not 
a local minimum solution but decays into the K* = 2~ 
solution after very long imaginary-time iterations (typi- 
cally 5 x 10 3 steps with At = 10 -3 MeV" 1 ). Reflecting 
different K v value of these two states, they are almost 
orthogonal to each other. To examine the physical re- 
ality of the K n = 0~ solution, we have performed the 
imaginary-time evolution of this solution with an extra 
constraint that the wave function is orthogonal to the 
K n = 2~ solution. In this procedure, we obtain a con- 
verged solution which is found to be almost the same as 
the original = 0~ quasi-stable solution without the 
orthogonalization. Because of these observations and a 
good correspondence to the measured spectra as shown 
below, we consider that the K v = 0~ solution is of phys- 
ical significance. 

Fig. shows density distribution of the internal Slater 
determinant $ in the three planes which include two prin- 
cipal inertia axes. The panels (a) are the SHF calcula- 
tion, and the panels (b) are the PPSHF for positive par- 
ity. They correspond to the ground state. The panels 
(c) and (d) are the PPSHF calculations for negative par- 
ity. The panel (c) is the lowest energy K w = 2~ state, 
and (d) is the K n = 0~. The nuclear shapes of (b) 
and (d) are characterized by the strong /3 30 deformation 
and are considered to reflect the a- 16 cluster structure. 
The clustering is stronger in the negative parity solution, 
which is consistent with the cluster model studies |(|. 
The lowest energy negative-parity solution shown in (c) 
looks prolate although it has a small j3%2 deformation. 

As for the lowest negative parity state with K n = 2~ , 
the excitation energy before AMP (the energy difference 
between the negative- and positive-parity solutions) is 
5.47 MeV and that after AMP (the energy difference be- 
tween = + and J 71 " = 2~ solutions in the parity and 
angular momentum projections) is 4.91 MeV. These val- 
ues are close to the measured 2~ excitation energy of 
20 Ne, 4.97 MeV. 

The next negative parity solution with K n — O - , 
which shows large fyo deformation, has a- le O cluster 
structure, as mentioned before. The excitation energy 
before AMP is 7.96 MeV and that after AMP is 6.42 
MeV, which are slightly higher than the measured ex- 
citation energy 5.78 MeV of the .V = \~ level. The 
different character of the two negative-parity solutions 
with = 2~ and K v = 0~ manifests clearly in the 
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FIG. 1: Density distribution of the internal state in Ne, in 
the xy, yz and zx planes where x, y and z axes represent the 
principal inertia axes, for (a) ground K 7 ' — + state of the 
SHF calculation, (b) ground K 71 = 0+ state of PPSHF, (c) 
excited K 77 = 2" state of PPSHF, and (d) excited K 77 = 0" 
states of PPSHF, respectively. The side of each panel is 14.4 
fm. The contour lines are plotted for every 0.02 fm -3 . 
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FIG. 2: Energy spectra of 20 Ne calculated with the AMP. 



spin-orbit energy, E^g, the expectation value of the two- 
body spin-orbit force. The K 77 = 2~ solution has large 
spin-orbit energy which is supposed to reflect dominant 
shell-model-like 5p-lh configuration. The spin-orbit en- 
ergy of K 71 = 0~ solution is much smaller, even smaller 
than that in the ground state. This is consistent with the 
development of a- 16 cluster structure in the negative- 
parity states which was seen in the density distribution 
in Fig.Hl 

In Fig. |3 we show the excitation spectra of 20 Ne af- 
ter AMP. For three rotational bands of K n = + , 2 _ , 



and - , calculated spectra are compared with measure- 
ments and with the results by the AMD (deformed AMD 
+ GCM calculation with Gogny force) Since each 

band is well characterized by the K quantum number, we 
simply show the energies given by E^^ — I^kk' ■ 
Although the band head energies of + , 2~, and are 
described reasonably well, the calculated moment of in- 
ertia deviates from the measured value. The calculated 
moment of inertia is too large for the ground state band. 
On the other hand, the calculated moment of inertia for 
the negative parity bands is slightly too small, opposite 
to the positive parity band. The bandhead energy of the 
K 77 = 2~ band in the AMD calculation is not as good 
as the PPSHF, probably because the present work has a 
better account of the single-particle wave functions. 

The AMD method gives better description for the high 
spin levels, 6 + - 8 + energy difference. In our calcula- 
tion, the AMP is performed from a single intrinsic state, 
while a change of nuclear shape is shown as the angular 
momentum increases in the AMD calculation [j^. This 
indicates that, if we incorporate the cranking in the PP- 
SHF framework, we might obtain better description for 
the higher angular momentum states. The discrepancy in 
the moment of inertia looks similar between PPSHF and 
AMD calculations. At present, we do not have a definite 
answer for the origin of this discrepancy. The pairing 
correlation ignored in the present calculation may be a 
possible answer. Bender et al. calculated ground state 
bands of some light nuclei in the HFBCS + GCM + 
AMP scheme Their calculation slightly underesti- 

mates the moment of inertia, opposite to our result. In 
the mean-field calculations, it has also been pointed out 
that the moment of inertia is sensitive to the time-odd 
component which we ignored 

In TABLE m we show observed and calculated intra- 
band B(E2) strengths of the K 77 = 0+ and K 77 = 2" 
bands of 20 Ne. For the sake of comparison, we show re- 
sults of the rotational model (Rot.) and the AMD + 
GCM 15]. Our calculation well reproduces B(E2) val- 
ues of the 2~ band. For the + band, our calculation 
somewhat underestimates the B(E2) values of low angu- 
lar momentum states. It should be noted that we do not 
introduce any effective charge. For 8 + — » 6 + transition, 
our result overestimates the B(E2) value, and is close to 
the rotational model. This is because we made AMP 
calculations from a single configuration. 

We next discuss 12 C. In the SHF calculations, one usu- 
ally obtain spherical ground state for a 12 C nucleus, al- 
though the rotational spectra is observed experimentally. 
The spin-orbit interaction favors spherical structure to 
gain energy through P3/2-P1/2 splitting. If one weak- 
ens the spin-orbit interaction slightly, oblate deforma- 
tion starts to appear in the ground state. We investigate 
change of the shape in the ground state by modifying the 
strength of the spin-orbit interaction, multiplying a con- 
stant factor xls to the two-body spin-orbit interaction 
of the Skyrme force. In Fig. |3 we show 0% and /?33 of 
the positive parity solution as a function of xls- The 0% 
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TABLE II: Observed and calculated intra-band B(E2) 
strengths of the K n = 0+ and K n = 2' bands of 20 Ne. 
For the sake of comparison, we show results of the rotational 
model (Rot.) and the AMD + GCM [p|. The observed val- 
ues are taken from Ref. [l8|. The values enclosed by curly 
brackets in the rotational model (rot.) are adjusted to exper- 
imental ones. 



K n = 


0+ 


B(E2) ft s 


Rot. 


AMD+GCM PPSHF 


2+^ 


0+ 


57 ± 8 


(57.0) 


70.3 


41.6 


4+ 


2+ 


71 ± 7 


81.4 


83.7 


59.9 


6+ -» 


4+ 


66 ± 8 


89.7 


52.7 


67.5 


8+ -> 


6+ 


24 ± 8 


93.9 


21.0 


75.1 


^ = 


2" 


B(E2)ob s 


Rot. 


AMD+GCM PPSHF 
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113 ± 29 


101 
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77 ± 16 
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FIG. 3: Deformation of C ground state as a function of 
xls which is a multiplicative factor for the two-body spin- 
orbit interaction. 



value in the ordinary SHF calculation is shown as well. 

We have found that, in the PPSHF calculation, an 
oblate shape with substantial deformation (triangle 
shape) appears if one employs a slightly weak spin-orbit 
interaction. This may be regarded as the appearance of 
the 3a clustering structure. The ^33 deformation does 
not appear in the ordinary SHF calculation. In addition, 
the oblate deformation starts to appear at larger xls 
value in the PPSHF calculation than in the SHF. In the 
ordinary SHF calculation, the appearance of the oblate 
deformation starts abruptly at about xls = 0.75 — 0.8. 
Note that, in order to remedy the isotope-shift problem in 
Pb isotopes, the spin-orbit potential should be weakened 
by about 30 percent from the ordinary one 0] . 

We have examined AMP for solutions with different 
xls values, and have found that the ground state rota- 
tional band of + , 2 + , 4 + is well described if we employ 
2:^5=0.8. At this value of x^g, we also made a calcu- 



TABLE III: Calculated results of 12 C at x LS = 0.8. The 
meaning of the listed quantities are the same as those in Ta- 
ble]!] 





A — U 


ncg£itiv6 p&rity 


qui? 


Energy (EXP) 


-92.162 


+9.64 


-92.162 


Energy (MF) 


-94.617 


+11.15 


-94.492 


Energy (MF+AMP) 


-96.499 


+9.69 




Els 


-13.229 


-7.550 


-15.191 




0.271 


0.638 


0.069 


1 


59.95° 


27.77° 


60.00° 


031 


0.0 


0.056 


0.0 


033 


0.327 


0.518 


0.0 
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FIG. 4: Density distribution of the internal state in C at 
xls = 0.8 for (a) SHF, (b) PPSHF (K w = + ), and (c) 
PPSHF (negative-parity) . See the caption of Fig. 0for expla- 
nation. 



lation for negative parity solution. We show below the 
results with this strength of the spin-orbit interaction. 
The results are summarized for energies and deforma- 
tions in Table II I II for intrinsic density distribution in 
Fig. 0] and for the energy spectra in Fig. [3] Calculated 
B(E2: 2+^0+) value is 5.14 e 2 fm 4 which should be com- 
pared with the observed value 7.8 ±0.4 e 2 fm 4 20]. 

As seen in Fig.|3](a), the nuclear shape is almost spher- 
ical in the ordinary SHF. At a slightly small xls value, 
oblate shape appears in the ordinary SHF calculation 
(Fig. In spite of the strong /? 2 and ^33 deformation, 
the spin-orbit energy is still large in the PPSHF calcu- 
lation, almost comparable to that in the spherical SHF 
calculation. This indicates that the closed p 3 / 2 configu- 
ration has still significant effects on the ground state. 

The negative parity solution shows strong /3a deforma- 
tion with triaxiality. It has also strong octupole defor- 
mation, mainly ^33 with small mixture of /3ai. Reflecting 
this mixture in shape, two configurations with different 
K quantum numbers also coexist in the solution. The 
panel (c) of Fig. 0] shows that the 3a clustering is much 



6 



> 

CD 



CD 
LU 



EXP CAL 



EXP 



CAL 



u it i- is* o- odd parity lowest 
r\ - 1 r\ - o energy state 



K"= + 



FIG. 5: Energy spectra of 12 C in the AMP. 



more developed in the negative parity than in the ground 
state. The AMP calculation gives reasonable description 
for excitation energy. Because of 7 deformation, not only 
3~ and AT states but also 1~ and 2~ states appear as the 
side band. 

It turns out that the structure of the negative par- 
ity solution in the PPSHF depends strongly on the xls 
value. If one employs the original Skyrme parameteriza- 
tion (xls = 1-0), we obtain K n = 1~ solution with P31 
deformation. On the other hand, for xls values smaller 
than 0.7, we obtain = 3~ solution with ^33 deforma- 
tion. At xls = 0.8, these two configurations mix up in 
the solution. Therefore, two distinct configurations co- 



exist in the negative parity at the excitation energies of 
about 10 MeV. In contrast to the 20 Ne case, two configu- 
rations do not separate but mix up. In such a situation, 
treatments beyond the present framework seem to be nec- 
essary, for example superposing multiple configurations 
in the generator coordinate treatment. We leave such an 
advanced treatment as a future problem. 

In summary, we propose the PPSHF method as a use- 
ful tool to study both excited and ground states simul- 
taneously. Self-consistent solutions of the excited states 
can be obtained in the negative parity, while the pos- 
itive parity solutions describe the ground state incor- 
porating certain correlation effects beyond the simple 
mean-field treatment. We show the feasibility of such 
calculations employing the uniform grid representation 
in the 3D Cartesian coordinate and achieving the three- 
dimensional angular momentum projection. The appli- 
cation to two N — Z nuclei, 20 Ne and 12 C, reveals that 
the obtained solutions show interesting deformations vio- 
lating reflection symmetries and incorporating clustering 
correlation. We will apply, for the future, the present 
framework for the systematic investigation and predic- 
tions of light nuclei, including exotic neutron rich nuclei. 

This work is supported by the Grant-in- Aid for Scien- 
tific Research (Nos. 14540369 and 14740146). A part of 
the numerical calculations is achieved on the supercom- 
puter at the Research Center for Nuclear Study (RCNP), 
Osaka University. 
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